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Neutrino speed experiments could be viewed not only as tests of Lorentz invariance but also as 
measurements of limiting propagation speed for all standard model species below certain depth 
where no direct metrological information is available. The latter option, hypothetically caused by 
some chameleon- type background, could be tested in the next installment of the neutrino speed 
experiments. We also show that that complementary constraints on the same class of models can be 
obtained with experiments testing clock universality in deep underground/underwater experiments. 
By considering the explicit QED model with particle-universal modification of propagation speed by 
a depth-dependent tensor background, we show that in general one should expect larger-than-GR 
shifts of the clock frequencies and clock non-universality. This can be tested by comparison of the 
narrow transitions in atomic clocks, which for a generic model could deliver a superior accuracy 
compared to the neutrino speed experiments. 

PACS numbers: 04.50.Kd 
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I. INTRODUCTION 



The OPERA experiment [l| created a stir of theoretical 
activity by suggesting that the speed of muon neutrino 
may be higher than the speed of light in vacuum, 



eOPERA 



^7, vac 



(2.37 ±0.32lg;|J) x 10" 5 , (1) 



with the significance of ~ 6a. These results agree with 
earlier measurement by the MINOS collaboration Q that 
also found e to be positive, consistent with (JT]), and yet 
consistent with zero at 2<r, 



cminos = (5.1 ±2.9) x 10" 



(2) 



While the claimed OPERA measurement has been shown 
to result from the initially undetected instrumental error, 
a new measurement by ICARUS Q is in perfect agree- 
ment with standard physics expectations, 



eicARUs = (0.01 ± 0.16 ± 0.37) x 10" 



(3) 



The initial results from OPERA led to the number of 
experiments testing the neutrino propagation speed, and 
more results are expected within the short time frame. 
Given that the neutrino speed experiments are reaching 
the 10 ppm (or better) accuracy of their measurements 
of c u , we ask the question whether such tests have addi- 
tional physics motivations, given that the original result 
(HJ) is incorrect. 

In this paper, we shall argue that the neutrino speed 
measurement provides an interesting probe of fundamen- 
tal constants at depths that are not yet accessible to di- 
rect metrological experiments. While we do not explicitly 
construct a dynamical model that could give e ~ 10~ 5 , 
we settle for an intermediate framework with the univer- 
sal density-dependent modification of the limiting prop- 
agation speed for all matter species. The purpose of this 



note is to point out that additional experiments with or- 
dinary atoms could also test deviations of c from c 7;V ac 
deep underground/underwater. 

Would continuation of the neutrino speed experiments 
provide nontrivial probes of Lorentz invariance? Despite 
last months numerous attempts to construct models that 
could "fit" the OPERA result, we note that none of the 
models with density-independent global modification of 
Lorentz invariance for neutrinos seem to work. Below 
we summarize the main arguments why large Lorentz 
violation (LV) for neutrinos is highly unlikely: 

1. As pointed out by Cohen and Glashow Q, the 
neutrino propagation faster than electron's limiting 
speed at the level of JT]) would result in rapid energy 
loss by neutrinos, which contradicts observations of 
energetic atmospheric and beam neutrinos. More- 
over, if speed of neutrinos is higher than the speed 
of quarks, a pion would not decay neutrinos above 
certain energies [5|. Note that the Cohen-Glashow 
argument applies not only to models that postulate 
global LV for neutrinos, but also to models that ex- 
ploit possible local modification of neutrino speed 
compared to other species caused e.g. by a tensor 
background Q . (For a divergent point of view see, 
however, Ref. 0.) 

2. Large global violation of Lorentz invariance for neu- 
trinos is in conflict with apparently normal timing 
of their arrival from the explosion of 1987a super- 
nova. 

3. Large global violation of Lorentz invariance for 
neutrinos is difficult to reconcile with tightly con- 
strained LV parameters for electrons, as the two 
sectors are connected by the W-boson loop [1, [§] . 

On account of these constraints, it is highly unlikely 
that Lorentz invariance is broken for neutrinos either 
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globally or locally at such a large level as e ~ 10~ 5 . 
Therefore, the direct measurement of c v does not pro- 
vide a superior test of Lorentz invariance for neutrinos. 
On the other hand, it is possible to speculate PUS!! that 
the limiting propagation speed is somehow modified for 
all species in the Earth's interior, and that neutrino speed 
experiments are the tests of this possibility. The neutri- 
nos in OPERA and ICARUS experiments propagate at 
average depth of ~ 6 km, reaching ~ 11 km in the middle. 
There are no direct tests of propagation speed performed 
at such depths, and therefore a possibility of the common 
to all particle change in c, however far-fetched, is not ex- 
cluded. Therefore, the results Q][3]) can be considered 
as tests of such possibility. One should also recall that 
gravitational field itself modifies the propagation speed 
for all species, but the modifications suggested by (UJ are 
many orders of magnitude larger. If indeed a propagation 
speed is modified at some depth below Earth's surface, 
one acquires two additional requirements for a model of 
this type: 

4. In-medium modification should sharply decrease 
near the Earth's surface in order to be in ac- 
cord with precise measurements of the gravitational 
force and the limiting propagation speed. 



5. Whatever backgrounds modify the propagation 
speed in the Earth's interior, they should not cou- 
ple to the 00 component of the stress energy for 
matter fields in order to avoid larger-than-gravity 
forces inside planets/stars. 

Models where in-medium properties of Lorentz invari- 
ant physical parameters such as masses and coupling con- 
stants are different from the same values in vacuum were 
introduced a few years ago [l2| ■ They are related to pre- 
vious ideas about the modification of the scalar-induced 
gravitational force by the presence of matter overdensi- 
ties [L3l4l5l ]. (Density dependence rather than redshift 
dependence could be an alternative interpretation of the 
non-zero result of Ref. [l6| that looked for the varia- 
tion of ckem in absorption systems at cosmological dis- 
tances.) Admittedly, models of in-medium modifications 
of the propagation speed are harder to construct, as they 
would require "condensation" of fields with non-trivial 
Lorentz indices. 

Assuming for a moment that some variants of density- 
dependence could modify the limiting propagation speed 
of neutrinos and all other species, we ask the ques- 
tion whether additional measurements performed with 
ordinary matter (not neutrinos) could test c and other 
constants deep underground/underwater. In the next 
section we discuss to what extent the (depth-induced) 
variation of the limiting propagation speed can lead 
to the variation of coupling constants and clock non- 
universality, and in the concluding section we propose 
new experiments that could detect such effects. 



II. LIMITING PROPAGATION SPEED AND 
CHANGING COUPLINGS 

The propagation speed of any matter species can 
be modified by the non-Lorentz invariant backgrounds. 
Consider a scalar field, with the Lagrangian modified by 
some tensor background h^: 



C <t> = 2^ 



1 

— m 

2 



V + h^d^dvcj). 



(4) 



Here we do not distinguish between upper and lower 
indices and perform the Lorentz summations with 
V^v — diag(l, — 1, — 1, —1). Moreover, we explicitly work 
in the system of units where h = c = 1, and the limit- 
ing velocity is 1 if the tensor background h^ v vanishes. 
The limiting velocity for the <f> particle travelling along z 
direction (E 3> m) is given by 



l-2h z 

1 + 2/ln 



1 - 2h 



00 



2h 2 



(5) 



One can see that negative /loo and/or h zz can "speed- 
up" the 0-field via effectively stretching time or short- 
ening distances. Besides the simplest possibility 
(|4]), one could introduce higher dimensional operators, 
h fl v\ K (j)di 1 d u d\d K 4> that will lead to the energy-dependent 
modification of the propagation speed. 

For the rest of this paper we use the following ansatz 



for the h^u field 



too 



hoi = 0; 



-e x diag(l, 1, 1), 



(6) 



so that (|4]) can be written as 
£<f = |(^t^) 2 — f(l + 2e)((9i</>) 2 , and we introduce 
the same modification for all fields of the standard 
model. The spatial anisotropy of e is not a crucial 
ingredient and is assumed for convenience. To comply 
with all requirements listed in the introduction, we 
take e = e(depth) to be some sharp function of the 
depth, with e = at the Earth's surface. The initial 
super-luminal claim by OPERA in our model would 
imply that at some critical depth zq one should assume 
that e deviates from and develops a positive value. 
The exact relation between e in (J6)) and eopera would 
depend on zq in a simple geometric way, but regardless 
of that measurement (fl]) would imply e > 2.5 x 10~ 5 , 
while of course ([3]) is perfectly consistent with e = 0. 

One could attempt replacing the tensor background 
(introduced here by hand) with some dynamical scalar, 
vector, or tensor fields: s, or H^ v . Going over 
to the canonical normalization of the kinetic terms 
for these fields, one can write down the interactions 
modifying the propagation speed of a generic SM field 
0, M' 4 (d^sd^) 2 ; M-^VM) 2 ; M^H^d^d^ !■ 
The immediate problem with the first two constructions 
is that in order to have any connection with OPERA re- 
sult, the scale M would have to be exceedingly low. The 
most recent model-building attempt to reconcile OPERA 
measurement with observations 11| employs a scalar 
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field, and M has to be below an MeV. Unfortunately, 
this is in plain contradiction with direct particle physics 
experiments. (E.g. electron-positron scattering remains 
consistent with the prediction of the SM to energies of 
^100 GeV, while the double-s exchange 1-loop diagram 
would lead to the significant modification of scattering 
for any M below O(10 GeV).) The tensor background of- 
fers perhaps the only reasonable hope for the dynamical 
model 6]. Still, even if we leave aside theoretical issues 
with UV completion, it appears difficult if not impossible 
to construct a weakly coupled model where dynamical 
H^ v (z) follows gravitational potential profile and does 
not run into contradiction with some observations. The 
model of Ref. @ that uses much enhanced coupling of 
H^ v to neutrinos faces a problem of Cerenkov radiation 
of electron-positron pairs, and an attempt to cure it by 
postulating the same interaction to electrons gives too 
much "anti-gravitational" force for electrons. For now, 
we shall assume that there is some consistent framework 
that leads to the in-medium condensation of along 
the lines of the proposals for the scalar field [HI, [I?} , al- 
though at this point it is an unproven assumption. 

Now we shall consider interacting fields and answer the 
question of whether the variation of e(z) could lead to the 
non-universality of clocks, or their abnormal speed-up or 
slow-down with depth, so that it could be picked up with 
dedicated experiments. (The connection between vary- 
ing c and coupling constants was previously discussed in 
Ref. [la].) For simplicity, let us consider the Lagrangian 
density of scalar quantum electrodynamics (QED) as a 
simplest model with gauge interactions, which we shall 
treat as a proxy to standard model. The unperturbed 
Lagrangian and its e-deformation are given by 

£qed - -if*, + +igA^\ 2 -™ 2 |$| 2 (7) 
Cc = V i-F^Fua + W» + igA^]*{d v + igA u )$) (8) 



i h l ^[6 1 (\D ll $\ 2 -m 2 m 2 )-8 2 -F? 



In these expressions, g is the gauge coupling and 
D fi = + igA^ is the covariant derivative, and trace 
is defined as h^^ = h^i]^. The two free parameters 5i 
and S 2 are meant to be order one, and they parametrize 
the model dependence. That is, in the first order in e 
their values do not affect the modification of the prop- 
agation speed. We now go to the ansatz ^ and take 
the adiabatic approximation where gradients of e are ne- 
glected. Then the sum of £qed and C e gives 

CQED+c = Ant + (|3q$| 2 - m 2 |$| 2 )(l + £&) 

-\d i ^\ 2 (l + e(2 + 6 1 )) (9) 
+ iE 2 (l + e(2- ( 5 2 ))-iB 2 (l + e (4-<5 2 )), 

where we separated terms bilinear in the fields from inter- 
actions. As evident from ©, the propagation speeds of 
photons and charged scalars are the same, c 7 = c$ = 1+e 
up to 0(e 2 ) corrections, and independent on ^1(2)- 



In order to determine whether one should expect ab- 
normal effects with clocks at O(e) level, we make redefini- 
tions of fields $ and A^ and distances dxi, while leaving 
the time variable unchanged. Selecting Aq = gauge, we 
have: 



(10) 



Making these changes in the action, S — J d 4 xC, and 
dropping primes over x, Aj, we read off a redefined 
equivalent Lagrangian to 0(e) level: 



$' = $ [1 + -(3 + Si)j ; dx' = dx(l - e) 
A' i = A i (l + ^-(5-S 2 



QED+e 



■^-m 2 |$| 2 



+\(d fl + ig(l-~(3-6 2 ))A fi )<Z>f 



(11) 



Thus, we have the same scalar QED theory, but the cou- 
pling constant is now changed to 



fl (l--(3-<S 2 ))j =a(l-ex (3-<5 2 )), (12) 



and is a function of depth, following the e(z) dependence. 
Notice that we do not have a change in mass of $ as a con- 
sequences of us choosing couplings of e to | c?o *& 2 — m 2 |$| 2 
combination. 

Different coupling constant means that the clocks build 
from "$-matter" working on the atomic transition of 
$ — $* bound state (~ a 2 m) will see the abnormal 0(e) 
difference when placed below zq. Also, different types of 
clocks with non-universal dependence on a will be sensi- 
tive to the e-induced change of frequencies. Elaborating 
on this, modification of couplings ([12")) implies two effects 
for the atomic clocks. Imagine a hydrogen atom that 
serves as atomic clock that uses optical and hypcrfinc 
transitions, ?ia; pt = |a 2 TO e c 2 and Tluim — 4 ^" e a 4 m e c 2 , 
where we used the proton magnetic g p factor, and have 
restored % and c for clarity. Then comparison of the 
optical or hyperfines frequencies on the surface and un- 
derground imply the following ratio of frequencies, 



u opt (z<zo) =1 _ 2e 
u opt {z = 0) 



o>hf(z < z ) _ ^ 
uitf(z = 0) 



■4ex (3-<y 2 ), 



(13) 
(14) 



while the comparison of optical and hyperfine transition 
by clocks underground should give a modified ratio: 



^hf 

Wopt 



(Z < Zq) 



32g p m c 
9m D 



a 2 (l - 2e x (3 - 6 2 )). (15) 



Using Eqs. p^|) - dT5l) as examples it is easy to predict 
what will happen for arbitrary clocks comparison since 
alpha dependence has been calculated for all known and 
proposed clocks -see e.g. 



Ref. 24 



Note that the de- 
pendence of the atomic unit of energy a 2 mc 2 on e (given 
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in Eq. (|13|) must be included when calculating depen- 
dence of the frequencies on e. This dependence was not 
included in all the previous calculations of alpha depen- 
dence in Ref. [24] since the results were presented in 
atomic units. For the purpose of the present work this 
dependence can be easily restored by adding factor a 2 
to all frequencies calculated in 24|- Then one should 
replace a by a e g from Eq. (fl"2"|) . 

Two things are further worth noticing: firstly, the 
modification of the limiting propagation speed for two 
species, and $, does not carry an unambiguous pre- 
diction for a e ff, as it depends on the free parameter 62 
not fixed by the requirement of the universality of c. On 
account of that the naive logic "c gets larger so that 
a = g 2 /(hc) gets smaller" does not hold. Secondly, there 
is a choice of £3 = 3 when the £qed = £qed+c- This 
choice corresponds to a situation when h^ v couples to the 
electromagnetic stress-energy tensor, h llv T llUl exactly as 
the linearized gravity would. In this case, to linear order 
in e, the clock universality will be preserved. 

Noting that e suggested by the original OPERA mea- 
surement |T]) is very large relative to the precision of mod- 
ern metrology, it is also interesting to investigate whether 
coupling of the background to the stress-energy operator 
would induce clock non- universality in 0(e 2 ) order. To 
that effect, we choose specific values of #1(2), 

£qed+£ = £qed + h^uTf^ 3, , (16) 

and follow the similar procedure to find that to 0(e 2 ) 
order this theory is equivalent to 

£ QED+e = ^(E 2 - B 2 (l - 4e 2 )) - m 2 |$| 2 

+ \(d»+*g(l-^)A,)<S>\ 2 . (17) 

Notice that the modification (fT6| creates 0(e 2 ) non- 
universality in the propagation speed of A^ and and 
this is why the e 2 -dependence persists for the free fields. 
Even if one neglects magnetic effects, the coupling con- 
stant is modified at 0(e 2 ) level. Only the full general- 
relativity-(GR)-like extension rj^ — > rj^ + of the 
original theory would preserve clock universality, which 
would entail additional (^(/i 2 ^) terms in (fT6|) with spe- 
cific coefficients. Thus, on the basis of (fTT|) and (|17|) . we 
conclude that barring a very special GR-like case, the 
modification of the limiting propagation speed for parti- 
cles leads to the clock non-universality. 



III. DISCUSSION: TESTING CLOCKS AT 
LARGE DEPTHS 

We have shown that one logical possibility - depth- 
dependent modification of c for all species - is not im- 
mediately ruled out by the variety of constraints on LV 
and by gravity tests. Moreover, such possiblity can be 
tested with the neutrino speed experiments. With the 



expansion of the baseline it is conceivable to probe the 
propagation speed as deep as 0(100 km). 

The purpose of this paper was not to built an explicit 
dynamical model with e.g. condensation of spin-2 fields, 
but to investigate whether depth-induced modification of 
the maximum propagation speed can be seen with more 
conventional means other than timing of neutrino events. 
We have shown that with a unique exception of pure 
GR-like coupling, one should expect an O(e) (or 0(e 2 ) 
in case of h liV T >lv coupling) deviations of couplings from 
their "surface" values, spatially linked to the deviation of 
propagation speed. An assumption that buidling a self- 
consistent dynamical models of this kind is possible is of 
course a "leap of faith", and several arguments can be 
presented why this is difficult to achieve. 

Over the years, there has been a concerted effort to 
test the GR theory in space [Hj]. Here we argue that 
the OPERA, MINOS and ICARUS results can be viewed 
as first metrological tests reaching 11 km underground. 
They give further insentives to test clock universality and 
their abnormal speed- up/slow-down at great depths. The 
GR-caused change in the frequency of clocks at depth h 
and on the surface is given by gh/c 2 and for h ~ 10 km 
is Auj/lli ~ 10~ 12 . Even the £oPERA" s ^ ze en?ec t is larger 
than GR shift by several orders of magnitude. Leaving 
aside the issue of the neutrino velocity, we believe that 
precision measurements underground are justified in their 
own right, as no systematic tests of this kind were ever 
performed. 

One could envisage two types of experimental set- 
ups to test the constants-changing-with-depth conjec- 
ture. First one could use the existing stable frequency 
emitters with known dependence on fundamental con- 
stants. Lowering them at great depths and comparing 
their frequencies either in-situ or by transmitting the 
signals to the surface will test clock universality. The 
second set-up uses two identical clocks synchronized on 
the surface, with one of them brought deep underwa- 
ter/underground for a period of time, with eventual com- 
parison of time measured by both clocks upon the return. 
Should any of such experiments indeed detect larger- 
than-GR effects of the depth on clocks, once could also 
experimentally determine zq, and further investigate pos- 
sible large gradient effects around z = zq. 

The deep underground locations could be used as a 
starting point for such tests. Indeed, deepest mines used 
for the underground science, such as e.g. Sudbury mine 
reach depths of 2 km. Similar depths in ice are reached by 
the IceCube collaborations operating at the South Pole. 
The deepest comercialy used mines in South Africa ex- 
tend 3.9 km underground. Ultimately, the "dream loca- 
tion" for such tests could be deep oceanic trenches and 
the deepest boreholes that extend as deep as 11 km or 
more (which is incidentally very close to the maximum 
depth along the OPERA and ICARUS neutrino trajec- 
tory). 

The connection between variations of the fundamental 
constants and the variation of the dimensionless ratios of 
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the transition frequencies of different atomic clocks is a 
well-researched subject. Current experimental sensitivity 
to the variation of a is better than 1 part in 10 16 |20l - l23l ] 
per year, far exceeding accuracy needs discussed in this 
paper. To check the link between the OPERA anomaly 
and variation of the fundamental constants it would be 
sufficient to compare commercially available Cs, Kb or 
quartz clocks which have accuracy 1CV 11 — 10~ 12 . The 
sensitivity of such clocks to the variation of the funda- 
mental constants has been calculated in j24j. Further 
six orders of magnitude improvement in sensitivity may 
be reached using optical clocks and the frequency comb, 
as well as with the use of atomic systems with closely 
degenerate levels (25j . 

Finally, we would like to comment that the density 



for many atomic clocks is very low, so that the ques- 
tion arises whether this would restore " surface values" for 
couplings inside atomic clocks deep underground. This 
question was answered in the models of chameleon- like 
varying constants [l2| , where it was shown that the real- 
istic model parameters ensure that the attainable sizes of 
the cavity are much less than the Compton wavelength 
of the chameleon field. In this case, the surface values for 
couplings will not be restored within the clock volume, if 
it is placed in an environment of non-zero e. 
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